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ABSTRACT 



in 



Aims. The availability of large amounts of multiwavelength data allows us to perform an extensive statistical analysis to look for corre- 
lations between different parameters of AGN. The physical implications of these correlations, when considered within the framework 
of current AGN models, can be enlightening in the resolution of the problems and the open issues which still characterise this class 

Oof sources. 
Methods. We presented CAIXA, a Catalogue of AGN in the XMM-Newton Archive, in a companion paper It consists of radio-quiet 
X-ray unobscured sources, which cover a range in X-ray luminosities between L2-iokcV = 2.0 x lO"*' and 3.9 x lO"*^ erg s"', and in 
redshift from z=0.002 to z=4.520. Here, a systematic search for correlations between the X-ray spectral properties and the multiwave- 
length data was performed for the sources in CAIXA. All the significant (> 99.9% confidence level) correlations are discussed along 
Q . with their physical implications on current models of AGN. 

5-H ' Results. Two main correlations are discussed in this paper: a) a very strong anti-correlation between the FWHM of the 11/3 optical 

' line and the ratio between the soft and the hard X-ray luminosity. Although similar anti-correlations between optical line width and 

, X-ray spectral steepness have already been discussed in the literature (see e.g., Laor et al. 1994, Boiler et al. 1996, Brandt et al. 1997), 

we consider the formulation we present in this paper is more fundamental, as it links model-independent quantities. Coupled with 
a strong anti-correlation between the V to hard X-ray flux ratio and the H/? FHWM, it supports scenarios for the origin of the soft 
excess in AGN, which require strong suppression of the hard X-ray emission; b) a strong (and expected) correlation between the X- 
ray luminosity and the black hole mass. Its slope, flatter than 1, is consistent with Eddington ratio-dependent bolometric corrections. 
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: such as that recently proposed by Vasudevan & FaMan (2009). Moreover, we cnt.ally review through vanous statistical tests the role 
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that distance biases play in the strong radio to X-ray luminosity correlation found in CAIXA and elsewhere; we conclude that only 
complete, unbiased samples (such as that recently published by Behar & Laor, 2008) should be used to draw observational constraints 
on the origin of radio emission in radio-quiet AGN. 

Key words. Galaxies: active - Galaxies: Seyfert - quasars: general - X-rays: general 

1 . Introduction increasing data archives. In particular, in the last few years large 

amounts of high-quality X-ray data have become easily accessi- 



' It is now well established that active galactic nuclei (AGN) are ^^j^ 
rS ] powered by accretion onto a central supermassive black hole 



. . powerea by accretion onto a central supermassive blaclc note , ^ . « ^ , «^xt ■ , 

$H St j — n 11 in/cn. r^ cc .a f j i We presented CAIXA, a Catalogue of AGN in the XMM- 

rrt ■ OLvnden- Bell 1969). Ditterent flavours ot unmcation models are » , . • • jtt 7-- . i ^„„„ 

'•^ , f ^ ■ , • • » f »u J u 1 Newton Archive, in a companion paper JBianchi et al.l2009 , B09 

quite successful in explaining most ot the varied phenomenol- ^ , . ' . ^ , ,. ^ : — 77 — T" 

f f J » 1 • u- u from now on). It consists ot all the radio-quiet X-ray unobscured 

ogy ot these sources in terms ot a fundamental scenario, which - '99 _9n » • ^ , • x7 , /.^xtx , 

, . , • . f A/^xT (Nh < 2 X 10 cm Active Galactic Nuclei (AGN) observed 

explains the characteristic signatures in the spectra ot AGN: a , . ,, . , , 

, ■• • .. , ■ u • 1 by XMM-A'ewfon in targeted observations, whose data are pub- 

thermal disc emitting in the optical and ultraviolet; a corona , • , • ■ , , 

cu ^ T ^ J • V u . • • . lie as ot March 2007. With its 156 sources, this is the largest 

ot hot electrons producing X-ray photons via inverse Compton , , ,-»^xrT 

„ . f u . f \u J- J . u- u catalogue ot high signal-to-noise X-ray spectra ot AGN. It cov- 

scattering ot the photons from the disc; dusty gas which re- ° . .? ° . . , , , ^^,^41 

, • . ■ f A ■ ■ A - ^ ers a range in X-ray luminosities between Li-iokeV = 2.0 X 10 

processes the nuclear emission into infrared emission and into , . „ ,„/tA 1 , . , , ^ ^ ^^r. 

V .1 . u . and 3.9 X 10^'' ergs"', and in redshift from z=0.002 to z=4.520. 

other X-ray spectral components. However, many aspects re- „ , • ^\ tx, « . „ . , 

^ ^ A ^ c A .1 c A/-^ivT ^ n * u f 11 For each source in CAIXA we have collected observable and 

lated to fundamental properties ot AGN are still to be fully un- ... 

A ^ u u A ^ A A physical parameters more directly related to the accretion onto 

derstood and many competing models have been developed and \^ ^ ,t,,,n 

, , , , r\ c -.f 1 .u- 1- J c- the supermassive black hole (BH), and less affected by contami- 

need to be tested. One fruitful approach to this kind ot inves- . , , , , r„, ^ , r-„T^ , 

^. . , , , J . ui nation by the host galaxy. Therefore, the EPIC pn X-ray data are 

tigation IS to gather as much multiwavelength data as possible , -' , , , . , , , ,. , • , , • 

, if f A e » f 1 complemented by multiwavelength data found in the literature; 

on large numbers ot sources and perform a systematic statistical „ „„,.,, ,^», • /T-,,„.»,^ ^,t^ 

, ° , , . . w u . A-ee . . BH masses. Full Width Half Maximum (FWHM) of Hfl, radio 

analysis looking tor correlations between ditterent parameters. ^^^^ tical fluxes 

Muc h progress have been do ne since the first seminal attempts ^ 

(e.g. lBoroson & Greenlll^ . flianks to the availability of ever- We refer the reader to B09 for aU the details on the selec- 

tion of the catalogue, the X-ray data reduction analysis and the 

Send offprint requests to: Stefano Bianchi coverage of the multiwavelength data. Here we investigate the 

e-mail: bianchiOf is .uniromaS . it correlations between the X-ray and the multiwavelength proper- 
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log(L,5VLj_,„) = {0.17±0.03) + (-0.98±0.09) log(H/Sp„„„ j„J [OLS B] 
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Fig. 1. Left: 0.5-2/2-10 keV luminosity ratio vs. H/3 FWHM. The 
anti-correlation between the two parameters is very significant. 
The two broken lines represent the 2000 km s ' boundary be- 
tween broad- and narrow-line objects and Lo.5-2/i'2-io = 1- See 
caption of Fig.|5]for details on the adopted symbols. 



ties of the sources in CAIXA. The physical implications of these 
results will be discussed within the framework of current AGN 
models, with particular interest on the open issues they still raise. 



2. Results 

We have performed a systematic correlation analysis between all 
the parameters which characterise the sources in CAIXA. The 
results are summarized in Table[T] In this paper, we consider sta- 
tistically significant only the correlations whose Null Hypothesis 
Probability (NHP) is lower than 1 x 10"^ (corresponding to a con- 
fidence level of 99.9%), with respect to the Spearman correlation 
coefficient p resulting from a censored fit (see Appendix lAl for 
details on this procedure, which takes into account both errors 
and upper limits on the Y variable). 

We will discuss in detail in the next sections all the signifi- 
cant correlations we found in CAIXA. 



2. 1. Correlations with 1-1/3 FWI-IM and the nature of the 
soft excess 

Looking at Table [T] the most significant correlation with the 
FWHM of Hj6 is that against the ratio between the soft (0.5-2 
keV) and the hard (2-10 keV) X-ray luminosity. The observed 
anti-correlation is shown in Fig.[T] 

Given the unknown origin of the soft excess in AGN, we in- 
troduced in B09 this X-ray luminosity ratio as a good parameter 
to characterize at least the basic property of the soft excess of the 
sources, i.e. its strength with respect to the primary X-ray emis- 
sion. As already noted in B09, it is very interesting to note that 
CAIXA does not contain narrow-line objects with a Lo.5-2/i2-io 



ratio lower than 1 and narrow- and broad-line objects are sig- 
nificantly different populations with respect to this parameter. 
However, just as no dichotomy appears in the optical properties 
between narrow- and broad-line objects, a smooth transition be- 
tween high and low X-ray luminosity ratios is also apparent from 
the data. The best fit for the anti-correlation takes this form: 



, ,^^0.5-2, 

log(- ) 

i-2-10 



(0.17 + 0.03) + (-0.98 + O.O9)log(Hj02ooo) (1) 



where ^82000 is the FWHM expressed in units of 2 000 km 
s~^. We remind the reader that this correlation, which involves 
the soft X-ray luminosity, was calculated on the CAIXA sub- 
catalogue which excludes objects with a measured column den- 
sity in excess of the Galactic one (see B09). 

The other correlation with H/3 found in CAIXA is the 
one with t he hard X-ray spectral slope, confirmin g previous 



claims (e.g. i Brandt et 



-ray spec tn 
al.l l 19971 If 



iconcelli et"ani2005 ). An anti- 



correlation between th e 0. 1-2.4 keV slop e and the li/3 line width 



i V slop e 

was also presented by Laor et al.' (1994") and later confirmed by 
iBoUer et al. (1996) and Laoretal. (1997). However, we note 
here that in CAIXA there is no significant anti-correlation be- 
tween the Hy6 and the soft X-ray photon index. This is likely an- 
other piece of evidence in favour of the fact that the values of F 
are very sensitive to the adopted models and it is therefore better 
to consider the correlations between more model-independent 
parameters. We therefore believe that our anti-correlation be- 
tween the X-ray luminosity ratio and Hp is a more basic physical 
correlation, because it links two a priori independent properties 
of radio-quiet AGN. 

The transition from broad- to narrow-line AGN, as well 
as the strength of the soft excess have been often interpreted 
as being due to a more fundamental physical paramete r, such 
as the accretion rate or the BH mass (e.g. Laor et al. 1997h . 
However, there is no significant correlation in CAIXA between 
Lo.s-ilLz-iQ and the Eddington ratio, the BH mass or any pa- 
rameter considered in this work other than the FWHM of H/3 (a 
part from the one with Fh discussed above). In other words, we 
do not have any indications that the strength of the soft excess 
may be related to any other important parameter, a part from the 
FWHM of H/3, not even the accretion rate. 

The 'classical' interpretation of the soft excess as direct 
emission from t he accretion dis c was questi oned when sev- 
eral studies (e.g. lGierlinski & Done 2004; Crumrny et al.||2006|) 
showed that the observed temperature of resulting black body 
is remarkably constant across orders of magnitude of luminosi- 
ties and BH masses. We confirmed the existence of a 'univer- 
sal' temperature, independent on luminosity or BH mass, on 
the large number of sources of CAIXA, if we substitute in the 
baseline model the soft X-ray power-law with a thermal compo- 
nent (B09). In the last few years, two main competing models 
have been proposed to explain the origin of the soft excess, both 
ascribing its apparent universal temperature to atomic physics 
processes. In one scenario, the soft excess arises from the en- 
hancement of reflection from the inner regions of the accre- 
tion disc due to light-bending effects, together with a strong 
suppression of the primary emission (e.g. Miniutti & Fabia^ 
2004). On the other hand, the soft excess can also be mimicked 
by absorption from a rel ativistically outflowing warm gas (e.g. 
iGierliriski & Donj|2004l). Regardless of th e details of the two 
models (and see also ' Merloni et al.l l2006b, for a self-consistent 
model of X-ray emission from inhomogeneous, clumpy discs), 
they both imply that the observed X-ray emission is suppressed 
with respect to the intrinsic one. 
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Table 1. Correlations between the main parameters which characterise the sources in CAIXA. 
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For each couple, the Spearman correlation coefficient p is reported (top), along with the logarithm of the NHP (bottom). If the latter is larger than 
-3.0, we consider the correlation with low significance and an 'X' is shown. Parameters with no significant correlations, namely EWgy, EWgciv 
and Fj are not listed in the table. Correlations involving Lo.5-2 or Lo.5-2/i-2-io were performed with the unobscured sub-catalogue of CAIXA. In 
this table, we also show trivial correlations for completeness, such as the ^-0.5-2/^-2-10 versus F;,, or the IIy6 FWHM versus Z-io//Z-£(/rf one. See B09 
and text for details. 



Therefore, the two models predict that the soft excess should 
be more prominent in X-ray weaker sources. As we have al- 
ready mentioned, there is no significant correlation between the 
X-ray luminosity ratio (our proxy for the strength of the soft 
excess) and any other parameter, apart from the FWHM of the 
HjS. We can only look for indications that the latter is correlated 
to any observed suppression of the intrinsic X-ray continuum. 
Indeed, we do find a significant correlation between the hard 
X-ray luminosity and the FWHM of HyS (see Table [T] an d Fig. 
|2]l, confirming previous claims (e.g. lPiconcelli et ani2005l) . And, 
more interestingly, we also find that the ratio between the opti- 
cal flux in the V band and the hard X-ray flux is significantly 
anti-correlated (p - -0.36 - NHP=2 x 10 with the width of 
the Hy6. This clearly shows that the correlation of L().5-2/i-2-io 
ratio with the H/3 FWHM is driven by the hard X-ray luminos- 
ity getting weaker, rather than the soft X-ray luminosity getting 
stronger. Therefore, there are observational indications that ob- 
jects with narrower H/? are X-ray weaker than normal broad-line 
sources. However, it is impossible from these results to disentan- 
gle between an absorption or reflection origin of the soft excess. 
Only detailed modeling of broadband spectra (incl uding high 
SNR data above 10 keV) can solve the issue (e.g. iPonti et all 
I2OO8.) . 

2.2. Correlations with Bi-i mass and tlie X-ray bolometric 
correction 

The strongest correlation with the BH mass in CAIXA is the one 
with the hard X-ray luminosity (see Table[T]and Fig.|3]l- This cor- 
relation is expected, because, if the accretion rate is a universal 
constant, the bolometric luminosity should rise linearly with the 
BH mass. Assuming that the bolometric correction, i.e. the ratio 
bet ween the bolome tric and the X-ray luminosity, is constant (as 
m lElvis et al.|[T994l) . flie correlation should be also linear with 



X-ray luminosity. The best fit slope for our data is slightly lower 
than 1, being 0.89 + 0.07, in agreement with previous results 
(Wang et al. 2009). This flatter slope could suggest that objects 
with larger BH masses have intrinsically lower accretion rates or 
are simply X-ray weaker with respect to their bolometric lumi- 
nosity. On the other hand, the correlation with the radio lumi- 
nosity ip - 0.59 - NHP=2 X 10"^ for the 20 cm catalogue) is 
consistent with a linear relation with the BH mass (see Fig. O, 
thus supporting the hypothesis that the total power of the AGN 
scales linearly with the BH mass, while high-luminosity objects 
may need a larger X-ray bolometric correction factor. 

The increasing optical-to-X-ray ratio with increasing opti- 
cal l uminosity is a well establi s hed result in AG N research (see 
e.g. lAvni & Tananbaumlll986t lJust etal.ll2007i and refer ences 
therein). This is in agreement with the X-ray bolometric cor- 
recti on factor incre a sing w ith luminosity adopted in CAIXA, 
after iMarconi et al.l (|2004|) . However, the correlation between 
the BH mass and the bolometric luminosity calculated in this 
way (p = 0.80 - NHP=7 x 10"^') gives a slope steeper than 
1, suggesting that the luminosity dependence of this correc- 
tion factor may be too stro ng at large luminosities. Recently, 
IVasudevan & FabianI (l2009l) have suggested that the X-ray bolo- 
metric correction may be dependent on the Eddington ratio, not 
the luminosity. We tried to adopt the average bolometric correc- 
tion reported in their paper for each Eddington ratio bin. The 
result is shown in the bottom right plot of Fig. [3] a correlation 
between the bolometric corrected X-ray luminosity and the BH 
mass with a slope consistent with 1 is recovered. 

Two other significant correlations with BH mass are found 
in CAIXA. The first is the correlation between the BH mass and 
z (see Fig. |4]i. This is expected: if the efficiency of the conver- 
sion of gravitational energy into radiative output is the same for 
aU AGN, the objects at high z, to be in CAIXA, must be lu- 
minous and, likely, with a large BH mass. Note that, at each 
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Fig. 2. Left panel: Hard X-ray luminosity vs FWHM ofH/3: the two parameters appear correlated. Right panel: Ratio between the 
optical flux (V filter) and the hard X-ray flux against the FWHM of Hy6: the observed anti-correlation suggests that narrow-line 
objects may be X-ray under-luminous. See text for details and caption of Fig.|5]for details on the adopted symbols. 



redshift, narrow-line objects have invariably the smallest BH 
masses: this is well known and some authors be lieve that the 
masses in these objects a re underestimated (e.g. iDecarli et alj 
120081; iMarconi et al.ll2008l) . This effect represents the other cor- 
relation with the BH mass found in CAIXA and can be explained 
in terms of the dependence of the BH mass on the H/? FWHM. 
This is of course due to the fact that most of the BH mass in 
CAIXA were primarily estimated from the H/S FWHM under 
the assumption of a virialised Broad Line Region (BLR). This 
would imply Mbh ^ FWHM(Hy8)^. The fact that the actual cor- 
relation in CAIXA requires a slightly steeper slope (2.96 + 0.23) 
is likely the result of the inhomogeneous provenience of H/3 data 
and methods for estimation of the BH masses. 

No other significant correlations with BH mass are found in 
CAIXA. In particular, the correlation with the Eddington ratio 
is lower than our threshold, even if larger than 99% confidence 
level (p = -0.29 - NHP=6 x 10"^). 

2.3. Correlations with radio emission and its origin in 
radio-quiet AGN 

In this section, we discuss the correlations between the radio 
properties of the objects in CAIXA and any other parameter. 
It should be borne in mind that the catalogue does not include 
radio-loud objects, as explained in detail in B09. 

In Fig. |5] (right panels), we plot the hard X-ray versus the 
radio luminosities (at 6 and 20 cm) for the objects in our cata- 
logue. The correlation is highly significant at both wavelengths 
(see Table[T]), with the following best fits: 

log(L6,4o) = (-0.48 + 0.08) -t- (1.19 + 0.08)log(L2-io,44) (2) 
log(L20,4o) = (-0.63 ± 0.06) + (1.13 + 0.06)log(L2-io.44) (3) 
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Fig. 4. BH mass vs redshift. See text for details and caption of 
Fig- E] for details on the adopted symbols. 



where radio and X-ray luminosities are given in units of 10"*" and 
lO'*"* erg s"', respectively. 

While the radio emission of radio-loud objects has been 
clearly estabhshed as originating in their relativistic jets (e.g. 
iBegelman et al.l[r984l) . the radio activity in radio-quiet sources 
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Fig. 3. Top 2-10 keV X-ray and 20 cm radio luminosity vs. BH mass for all the sources in the catalogue with an estimate of the 
latter. A linear regression fit is superimposed on the plots. In the left panel, the curves referring to Eddington rates of 1 and 0.01 
£ire also shown. Bottom Bol ometric luminosity (left: luminosity-dependent ( Marconi et al.ll2004l) : right: Eddington ratio-dependent 
dVasudevan & Fabianll2009l) ) vs. BH mass for all the sources in the catalogue with an estimate of the latter A linear regression fit is 
superimposed on the plots. See text for details. 



is still puzzling. Several authors have shown that X-ray 
and radio luminosity stron gly correlate in radio-quiet AGN 
(e.g . Brinkmann et al. 2000; iPanessa et al. l l2007HLiet al1l2008t 
iLaor & Behar.,200 8). This result is generally interpreted as ev- 
idence that also in radio-quiet objects radio emission can be 
considered a direct tracer of the nuclear activity, either through 
a common accretion process or through the scale-invariant de- 



pen dence of th e radio-emitti ng jet on t he accretion physics (e.g. 
.Heinz & Sunvaev 2003; .Merloni et al.ii,200 3 ) . 

The presence of the strong correlation found in CAIXA is 
in agreement with previous results, even if the slope of our best 
fits (around 1.1-1.2) is somewhat steeper than the linear rela- 
tion s found by Br inkmann et al. (2000), Panessa et al. (200^ 
and iLaor & Behad ( i2008h for radio-quiet objects. Our slope re- 
mains steeper even if we consider the correlation between the ra- 
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Fig. 5. Radio to X-rays plots. From top to bottom: 2-10 keV vs 6 cm emission, 2-10 keV vs 20 cm emission. Flux-flux correlations 
are on the left, luminosity-luminosity correlations on the right. In each case, the analytical expression for the best fit is on the 
top. The different symbols refers to the classification of the objects, on the basis of their absolute magnitude and H/3 FWHM: 
NLSY, narrow-line Seyfert 1; BLSY, broad-line Seyfert 1; NCSY, not-classified Seyfert 1 (no Hj8 FWHM measure available); NLQ, 
narrow-line quasar; BLQ, broad-line quasar; NCQ, not-classified quasar (no H/3 FWHM measure available). 



dio emission and the soft X-ray luminosity (see Table[T|and Fig. 
6t, as the above rn entioned authors do (with the exception of 
Panessa et a n i2007h . Other authors found correlations far from 
being linear, and significantly different between radio-loud and 
radio-quiet objects (e.g. Li et al. 2008). We therefore decided to 
investigate possible biases affecting this correlation. 

It is known that care must be taken in assessing the real sig- 
nificance of any luminosity-luminosity correlation and in giving 



a physical meaning to the best fit slope, because of the effects in- 
troduced by distance and the selection of the catalogue. Indeed, 
as expected, the luminosity is tightly correlated to the redshift, 
due to the sensitivity limit of our observations (although it varies 
depending on the exposure time of each observation). 

In the left panels of Fig. |5j we show the same correlations 
as above, but with fluxes instead of luminosity. The correla- 
tions are still statistically significant, but the significance is lower 
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Fig. 6. Radio 20 cm vs. soft X-ray luminosity, along with best 
fit. See caption of Fig.|5]for details on the adopted symbols. 



{p = 0.49 and 0.67, with NHP=7 x 10^^ and 5 x lO'^^ for Fe 
and F20, respectively) and the scatter is much larger than that be- 
tween luminosities. Moreover, the best slopes are flatter than the 
ones found with luminosities, in particular for the measures at 
20 cm, where there are many upper limits. This suggests that the 
role of distance in the previous luminosity-luminosity plots is 
indeed important. The reason is obvious: only objects with high 
luminosities in both bands can be seen at large distances, while 
both high-and-low luminosity objects can be seen at small dis- 
tances (but the small volume often reduces the numb er of hig h- 
L ob jects). This was shown in the simulations by llsobe et al.l 
(11986) . who demonstrated that, if the upper limits are kept, then 
survival analysis methods can in principle recover the correct re- 
lationship. However, our catalogue can still be affected by this 
effect, because, even if we keep upper limits (both on radio and 
X-ray fluxes), neither the X-ray nor the radio data-sets are by 
any means complete and they are selected only on the basis of 
the SNR of the first, thus likely introducing a bias towards high 
flux sources. We therefore do not consider the existence of a 
significant correlation between fluxes a proof significant enough 
that the distance bias does not play any role in the luminosity- 
luminosity correlation in CAIXA. 

Accordingly, we performed a partial Kendall t correlation 
test on censo red data, using the distance as the third variable, as 
proposed by Akritas & SiebertI (Il996i) . Also taking into account 
the correlation with distance, the X-ray versus radio luminosity 
correlation still remains significant in all the cases: coefficient 
T = 0.24 and variance cr = 0.056 (NHP=3 x 10"^) for 6 cm, 
T = 0.22 and a = 0.046 (NHP=1 x 10"'') for 20 cm. 

We finally performed the "scrambling test" (e.g. iBregmarJ 
l2005HMerloni et al.l2006al) . The principle at the heart of this test 
is very simple: we keep together each couple of X-ray luminosity 
and redshift and "scramble" the radio fluxes, assigning each one 
to a random L2-10/Z couple. Then, adopting the newly assigned 
redshift, we calculate the radio luminosity from the random flux 
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log(z) tog(z) 

Fig. 9. Redshift distributions of the sources in CAIXA as a func- 
tion of hard X-ray flux. The distributions clearly evolve with 
flux, revealing an intrinsic bias of the catalogue, due to selec- 
tion effects. See text for details. 



and correlate it to the assigned X-ray luminosity. Any physical 
correlation intrinsic to the original data-set must disappear in the 
new artificial data-sets, because the radio flux of a given object 
is now associated to the X-ray luminosity of another one. We 
produced 100000 scrambled data-sets with this procedure for 
the 6 and the 20 cm radio fluxes and we calculated for each of 
them the Spearman correlation coefficient, the Kendall partial 
correlation coefficient and the censored fit. The results are shown 
inFig.|7] 

It is striking to see how the Spearman correlation coefficients 
for the real data-sets are well within the distribution of the simu- 
lated ones: in the case of the 6 cm data, around 87 000 values of 
p are actually larger than the 'real one' . Therefore, this test sug- 
gests that the probability that the original correlation is driven 
entirely by distance is very high, i.e. 87%. This probability is 
much lower for the 20 cm data (which includes upper limits), but 
still large: 8%. The simulated distributions of the Kendall partial 
correlation coefficients are instead shifted towards less signifi- 
cant values with respect to the real ones, but still the probability 
that distance plays the main role is rather high, being 7 and 4%, 
for 6 and 20 cm, respectively. However, in Fig. [8]we also show 
the distributions of the best fit slopes for the simulated data-sets 
compared to the "real" ones (right plots): the latter are much 
flatter. This indicates that the observed slopes are not random, 
and may still reflect a real intrinsic correlation. On the other 
hand, this effect may also be due to the fact that CAIXA con- 
tains few objects with large redshifts and no radio fluxes below 
a certain limit. These two factors contribute to enhance the pop- 
ulation of very luminous radio objects in the simulated data-sets 
(when an average or large radio flux is associated to a large z), 
while the region with high redshift and low radio luminosity is 
still under-populated, leading to steeper slopes, i.e. objects more 
radio-luminous than the real ones. 
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Fig. 7. Scrambling test results for the correlations between the hard X-ray luminosity and the radio luminosity at 6 cm (upper plots) 
and 20 cm (lower plots). The left panels show the distribution of the Spearman coefficients for the 100000 simulated data-sets, in 
comparison to the one measured for the real data-set (dashed vertical line). Same in the right panels, but for the Kendall partial 
correlation coefficient t divided by the variance cr. 



As a last check, we also performed the "scrambling test" for 
the flux-flux correlations. In this case, the simulated Spearman 
correlation coefficients are centered (both for 6 and 20 cm) in 
0, as expected for random uncorrected points (see FigfTOli. Not 
even a simulated data-set has a correlation coefficient as large as 
the ones measured in the real data-sets, in agreement with the 
very low estimated NHP in the latter However, the strong corre- 
lation found between the X-ray and radio fluxes in CAIXA may 
still be implicitly due by distance effects. Due to selection ef- 



fects. X-ray fluxes are tightly linked to redshift. This is clearly 
shown in Fig. |9] where we plot the redshift distributions of the 
sources in CAIXA in four hard X-ray flux bins. The evolution of 
the distributions is evident: faint objects are clustered at high z, 
while bright sources are peaked at low z. Therefore, although we 
cannot exclude the presence of a physical correlation between 
hard X-ray and radio fluxes in CAIXA, we are unable to disen- 
tangle it from the selection effects present in the catalogue. 
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Fig. 8. Scrambling test results for the correlations between the hard X-ray luminosity and the radio luminosity at 6 cm (left) and 20 
cm (right). The best fit slope distributions for the 100000 simulated data-sets are shown, in comparison to the ones measured for 
the real data-set (dashed vertical line). See text for details. 



In conclusions, the observed correlation between radio and 
hard X-ray luminosity appears extremely significant on the basis 
of standard statistical tests (Spearman correlation coefficients, 
partial Kendall correlation), but the role of distance effect may 
be dominant in determining the significance and functional rela- 
tion. In particular, the scrambling test shows that even the widely 
used partial Kendall correlation may underestimate the distance 
effects, thus care must be taken when interpreting the results de- 
rived from it. However, random scrambling cannot reproduce the 
observed slope of the luminosity correlation, and destroys the 
flux-flux correlation, which is instead observed to be significant. 
Both results may represent the revealing signs of a real underly- 
ing correlation, but they can also be induced by the incomplete- 
ness of CAIXA. 

Therefore, although a physical correlation may be present in 
our data, we conservatively cannot draw any conclusion from 
the X-ray/radio correlations in CAIXA. We would like to stress 
that this is true for CAIXA, which includes upper limits, but it 
is by no means a flux-limited or volume-limited sample, while 
distance effects may of course play less important roles in other 
catalogues or in complete samples. 

In the left panel of Fig. [TT] we plot the correlation between 
the two radio-loudness parameters used to exclude the radio- 
loud objects in CAIXA, one based on the optical fluxes found 
in literature, the other on the X-ray fluxes measured in this work 
(see B09 for details). As expected, a highly significant linear re- 
gression between the two parameters is found (see Table[Tli. The 
slope of the best fit (0.88 ± 0.06) is p erfectly in agreeme nt with 
the on es found bv iTerashima & Wilson ( 2003) and Panessa et al.l 
(l2007h . even if their samples include also radio-loud objects, 
thus spanning a larger range both for Rx and Rq. 

In the right panel of Fig. (TT] we show the X-ray radio- 
loudness parameter versus the Eddington ratio. Ho (2002) found 
an highly significant anti-correlation between these two pa- 



rameters, suggesting that weakly active nuclei are powered by 
advection-dominated accretion flows (ADAFs), which give rise 
to radio-loud spectra. Although this trend seems to be observed 
in CAIXA, the anti-correlation is not significant (p - -0.14 - 
NHP=0.19). However, this is not unexpected, because our cata- 
logue does not include radio-loud objects and, indeed, the large 
Rx-low Lboi/Ledd region is not populated in our plot. Finally, we 
note here that no correlation is apparent between the soft to hard 
X-ray luminosity ratio and the radio luminosity. 

2.4. The X-ray spectral index 

Apart from the correlations with Hfi and the X-ray luminosity 
ratio, already discussed in Sect. 12.11 there are no other signifi- 
cant correlations with the hard X-ray spectral index in CAIXA. 
In particular, in the literature there are claims of a correlation 
between F/, and the 2-10 keV luminosity (e.g. iDai et al]|2004t 
Saez et al. 2008), but there are also se veral studies that do not 
find s uch a correlation (e .g. lReeves & Turner 2000.;, George et al.l 
20001: 1 Winter et al.ll2009l) . In CAIXA, fliis correlation is not sig- 
nificant (96% confidence level): we note her e that this confi- 
den ce level is not far from the ones quoted bv lDai et~aLl (|2004|) 
and Saez et alj (1200 8). unless the correlation is calculated in par- 
ticular bins of (relatively high) redshift. In other words, the low 
significance of the F/i/La-io correlation in CAIXA is in agree- 
ment with previous studies, even with those which claim the 
correlation is real. Nonetheless, we remind that significant dif- 
ference between the spectral index of Seyferts and quasars is 
measured in CAIXA (Tsy = 1.66 + 0.05: Fqso = 1-80 + 0.05: 
B09). We showed in B09 that the introduction of a luminosity- 
dependent Compton reflection component to the baseline model 
does not remove the difference between these two populations. 
Further tests should be performed to understand the real driver 
for this difference. 
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Fig. 10. Scrambling test results for the correlations between the hard X-ray flux and the radio flux at 6 cm (left) and 20 cm (right). 
The distributions of the Spearman coefficients for the 100000 simulated data-sets are shown, in comparison to the ones measured 
for the real data-set (dashed vertical line). 
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Fig. 11. Left: Radio-loudness parameter vs. X-ray radio-loudness parameter (see text for definitions) for all the sources in the 
catalogue with a radio measure. The dotted horizontal line represents the canonical separation between radio-loud and radio-quiet 
quasa rs, while the dashed lines the separation for Seyferts, based both on the optical as well as on the X-ray radio-loudness parameter 
(after iPanessa et al ]l2007l) . A linear regression fit is superimposed on the plot, along with the best fit parameters. We refer the reader 
to B09 for the detailed selection criteria adopted in CAIXA with respect to the radio-loudness. Right: X-ray radio-loudness parameter 
against Eddington ratio. Although a trend is observed, the anti-correlation is not statistically significant. See caption of Fig. |5] for 
details on the adopted symbols. 
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Similarly, we do not find a significant correlation between 
the photon index and the BH mass or the Eddington rate, but the 
significances found by Piconcelli et al. (2005) are lower than our 
threshold, and the results recently found by iKellv et al.l (l2008l) 
are strongly dependent on how the BH masses are calculated. 
O n the other hand, the cor relation with the Eddington ratio found 
bv lShemmeret all ( 1200 8') is more significant than our threshold, 
but only when high-z and high-luminosity quasars are taken into 
account. 

2.5. The iron lines 

iBianchi et al.l (|2007|) already discussed in detail the strong anti- 
correlations between the Equivalent Width (EW) of the neu- 
tral iron narrow emission line and the 2-10 keV luminosity and 
Eddington rate. Despite the exclusion of a source (PMNJ0623- 
6436) in CAIXA, the results published in Bianchi et al.1 (l2007h 
do not suffer any variation. Recently, .Dadina C200S ) confirmed 
the so-called 'Iwasawa-Taniguchi' (IT) effect considering the 
20-100 keV luminosity, with a slope fully consistent with the 
one found with CAIXA in the 2-10 keV band. 

It is interesting to note that the likely interpretation of the IT 
effect as due to the covering factor of the torus decreasing with 
luminosity would imply a decrease of the fraction of Compton- 
thick objects with luminosity. This is indeed what found by 
iDella Ceca et al. (2008). Once extracted from their results the 
fraction of Compton-thick objects with respect to the total AGN 
population, which can be considered a proxy for the covering 
factor of the torus, we found that it decreases with luminosity 
following a slope of ^ -0.22, in very good agreement with the 
IT effect. 

We note here that the IT effect is also in agreement with 
the recent discovery of a non-linear relation between thermal 
emission from dust and optical luminosity in AGN, which im- 
plies a decrease of the cov ering factor of dust with luminosity, 
with a slope of ^ -0.18 ( Maiolino et al.ll2007i) . i.e. very close 
to the one which characterises the IT effect. Although the agree- 
ment is very good and deserves further investigation, it must be 
noted that the two methods sample, in principle, different ma- 
terials. The iron Ka line is mainly produced in Compton-thick 
gas, while the dust can be associated both to Compton-thin and 
thick gas. 

As a final note, we do not find any significant correlation 
involving the EWs of the ionised iron hnes. 

3. Conclusions 

In this paper, we investigated the correlations between the X-ray 
and the multiwavelength properties of the sources in CAIXA, 
a Catalogue of AGN in the XMM-Newton Archive, presented 
in detail in a companion paper (B09). The main results can be 
summarised as follows: 

- Correlations with H/3 FWHM. A very strong anti-correlation 
between the FWHM of H/3 and the ratio between the soft 
and the hard X-ray luminosity is present in CAIXA. In par- 
ticular, our catalogue does not contain narrow-line objects 
with a Lo.5-2/i2-io ratio lower than 1. Although this correla- 
tion reflects the same physical driver as correlations already 
published in th e literature between the soft X-ray slope and 
the H/3 FWHM dBoUer et alJl996tlLaor et aiTT99 7). we con- 
sider that it bears a more fundamental meaning, because it 
links two more model-independent quantities (it is interest- 
ing to notice that no correlation between the soft X-ray slope 



and the H/3 FWHM is found in CAIXA). The ratio between 
the optical flux in the V band and the hard X-ray flux is also 
significantly anti-correlated with the Hj6 line width, suggest- 
ing that objects with narrower H^S are X-ray weaker than nor- 
mal broad-line sources. 

- Correlations witli the BH mass. The strong correlation be- 
tween the X-ray luminosity and the BH mass has a slope 
flatter than 1, suggesting that high-luminosity objects may 
be X-ray weaker A luminosity dependent bolometric fac- 
tor could exp l ain th is result, but the one proposed by 
iMarconi et al.l d 2004 seems too strong at high luminosi- 
ties. On the other hand , with the bolometric factor pro- 
posed by V asudevan & Fabian (2009), which depends on the 
Eddington ratio, a relation with a slope consistent with 1 is 
recovered. In any case, that the intrinsic total power of AGN 
indeed scales linearly with BH mass is suggested by the lin- 
ear relation between the radio luminosity and the BH mass. 

- Correlations with radio emission. In CAIXA there is a very 
strong correlation between the radio (6 or 20 cm) and the 
X-ray (soft or hard) luminosity. Similar correlations were 
presented and discusse d by several authors in the past (e.g . 
iBrinkmann et al J 120001; IFanessa et al.1 120071: iLi et"ani2008l) . 
We have critically reviewed the effect that the incomplete- 
ness of the sample may have on this correlation. The "scram- 
bling test" suggests that indeed a distance bias may be domi- 
nant, even if the partial Kendall t test often used in the litera- 
ture to account for this bias suggests otherwise. Our conclu- 
sion points to the need of using co mplete unbiased samp les, 
such as that recently discussed by iLaor & BehaJ (l2008l) . to 
draw observational constraints on the nature of radio emis- 
sion in radio-quiet AGN. 
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Appendix A: Linear fits on censored data 

This appendix is very similar to the one presented in IBianchi et alj j2007l) . We 
include it here for the reader's convenience and to report the inclusion of the 
scattering en'or in the overall uncertainty estimated for the fit parameters. 

Since we are dea ling with a large number of upper limits, we followed 
IGuainazzi et all )200d) to perform a linear fit, taking into account both upper 
limits and errors on the y pai'ameter Their m ethod is an exte nsio n of the re- 
gression method on censored data described bv lSchmitH (1985") and'Is obe et alj 
( 1986) and can be summarized as follows. Several Ordinary Least Squares Y 
vs. X (OLS Y/X) or Bisector (OLS B) fits (each plot reports on the upper right 
which method was used) was performed on a set of Monte-Carlo simulated data 
derived from the experimental points according to the following rules: a) each 
detection was substituted by a random Gaussian distribution, whose mean is the 
best-fit measurement and whose standard deviation is its statistical uncertainty; 
b) each upper limit U was substituted by a random uniform distribution in the 
interval [0,U]. The mean of the slopes and the intercepts derived from the fits of 
each data set are our 'best fit' . Note that errors and upper limits on the x param- 
eter are not considered in the procedure, although plotted. This does not affect 
the results significantly, given the fact that they are typically much smaller than 
the errors on the y. Moreover the errors on the best fit parameters only take into 
account the experimental uncertainties, being the dispersion of these parameters 
within the Monte Carlo simulated data sets. Any intrinsic, physical scattering of 
the data points is not included. To test the significance of this anti-correlation, 
we calculated the Spearman's rank coefficient for each Monte Carlo simulated 
data set: the mean value is taken as the Spearman's rank coefficient of our 'best 
fif . 
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